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Can hidden correlations mimic a variable fine structure constant? 
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Abstract. Murphy et al. (2003a, MNRAS, 345, 609) claim to find evidence of cosmological variations of the fine structure 
constant a in the spectra of intervening QSO absorption line systems. We find that this result is affected by systematic effects. 
The a values estimated in individual line systems depend on the set of atomic transitions used and therefore the quoted de- 
pendence on the cosmic age may reflect the fact that different sets of transitions are used at different redshifts. A correlation 
between line shifts and relative optical depths of the atomic transitions is also present. This correlation is very tight for a high- 
redshift subsample and accounts for the anomalous dispersion of the a values found by Murphy et al. (2003a) in this subsample. 
The above correlations are consistent with a scenario in which gravitational redshift, caused by the gravity of the dark halo, 
contributes to the shift in frequency of individual components. Gravitational redshift causes differential spectral shifts of the 
same order as magnitude of those measured. In the presence of line misidentification, these shifts can be interpreted in terms of 
a variable a. In order to verify the gravitational redshift hypothesis, a direct access to Murphy et al. (2003a) data, or to a large 
amount of new high resolution data, is necessary. 
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1. Introduction 



The hypothesis of a cosmological evolution of the fine struc- 



ture constant 


a has been investi 


gated in 


a series of pa- 


pers (Murphv, Webb & Flambaum, 


2003a; 


Webb et al., 


2003; 


Murphv et al. 


l2001aUWebb et alll200 ill 19991 also collectively 



referred to as MWF), by applying a sophisticated spectral anal- 
ysis, called the "many multiplet" method (hereafter MM), to 
large samples of absorption line systems in QSO spectra. This 
method is based on a^ 2 analysis of multiple-component profile 
fits in several transitions, and relies on theoretical calculations 
of how the relative positions of different transitions change due 
to a variations (e.g. lDzuba et alll2002i) . The use of many tran- 
sitions improves the accuracy of a determination in each ab- 
sorption system, compared to fits of only two transit ions (like 
the alkali-doublet method, e.g. Murnh vet all 120013) . In addi- 
tion, the use of large samples entails a much smaller quoted un- 
certainty on the average v alue of Aa/a ( up to 1 28 absorption 
systems were analyzed in Murp hy et all l2003a). The conclu- 
sion of MWF is that the value a was slightly smaller at earlier 

cosmological epochs (Aa/ a 0.5 X 10~ 5 , with a 5-cr s ignif- 

icance). Ancillary papers llMurphv et all l2001bl |2003b) have 
reached negative conclusions on potential systematic effects 



which could mimic a variable fine structure constant (hereafter 
VFSC). 

A recent work bv lChand et all d2004l) . based on more ac- 
curate measurements in high quality data relative to a smaller 
sample objects, uses the MM method with well-defined selec- 
tion criteria. No cosmological variation of Aa/a is found, with 
a 1-cr uncertainty of 0.06 x 10~ 5 (namely an order of magni- 
tude below the values quoted by MWF). Chandetal 
avoid very weak and blended lines and point out that when one 
does not do so, spurious detections are frequently seen. U sing 
a different but very robust method. fBahcall et all (|?004) also 
do not fi nd any relevant time dependence of the fine structure 
constant. Bah call et all 12004) point out several inconsistencies 
in the results of MWF which imply that systematic uncertain- 
ties due to misidentification of lines might be significant in the 
MM method used by MWF. In order to directly test the hy- 
pothesis of component misidentifications, one needs to investi- 
gate in detail the original data and the routines used by MWF 
in their analysis: this information is, how ever, not acces sible 
to the community. As e mphasized bv lChand et all ll2004l) and 
bv lBahcalletall ll2004 . MWF have not described their algo- 
rithm of line identification, the confidence level adopted and the 
physical assumptions made. So, line misidentifications cannot 
be excluded in their results. 



Send offprint requests to: R. Bandiera, 
e-mail: bandiera@arcetri . astro . it 



Given the overwhelming consequences in the fundamen- 
tal physics and cosmology of the variability of the fine struc- 
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ture constant (see e.g. iMartinsl 120031 lUzanl 120031 12004 . any 
newly suggested potential source of bias in the determination of 
Aa/a must be examined in detail. Direct access to MWF data 
and analysis would facilitate this process. At the moment, how- 
ever, we are forced to approach the problem in different ways. 
In the present paper we investigate t he stati stical properties of 
the data presented by Murph y et al.l J2003al their Table 3), and 
we find correlations that can be explained only in terms of a 
systematic effect, not accounted for by MWF (Section 2). In 
Section 3, we put forward the hypothesis that differential grav- 
itational redshift (hereafter GRS), within absorption systems, 
may be responsible for the correlations found and of the non- 
zero value of Aa/a. An important feature of GRS, with respect 
to other effects, is that symmetry between red- and blueshifts is 
broken, and therefore a systematic effect on the average values 
of the spectral shifts may be accounted for. This effect is very 
tiny, and it can be detected only under particular conditions. 
In the MM analysis, GRS appears in only if misidentifications 
of some components are present. In Section 4, we search for 
correlations between a and some indicator of the gravitational 
potential, for various absorption systems: the presence of such 
correlations represents a strong clue that some differential red- 
shift is of gravitational origin. Our last Section summarizes and 
concludes that GRS is a possible candidate for the origin of the 
measured differential redshift. 



Table 1. Atomic data. 



Ion 


log Mel 


Wavelength 


fx 10 3 


logr 


Qx 10 2 


Mgi 


-1.6520 


2852.96 A 


1810.0 


-0.8213 


0.2454 


Mg ii 


0.3372 


2796.35 A 


612.3 


0.6885 


0.5900 






2803.53 A 


305.4 


0.3876 


0.3364 


Alii 


-1.0454 


1670.79 A 


1880.0 


-0.4306 


0.4511 


Aim 


-1.5818 


1854.72 A 


539.0 


-1.4641 


0.8606 






1862.79 A 


268.0 


-1.7657 


0.4024 


Sin 


0.4659 


1526.71 A 


127.0 


-0.1288 


0.1038 






1808.01 A 


2.18 


-1.8207 


0.9601 


Cm 


-1.6784 


2056.26 A 


105.0 


-2.2263 


-2.2763 






2062.24 A 


78.0 


-2.3542 


-2.5799 






2066.16 A 


51.5 


-2.5336 


-2.7563 


Fen 


0.0000 


1608.45 A 


58.0 


-0.9124 


-1.9301 






1611.20 A 


1.36 


-2.5415 


1.6918 






2344.21 A 


114.0 


-0.4553 


2.9396 






2374.46 A 


31.3 


-1.0111 


3.8941 






2382.76 A 


320.0 


0.0000 


3.5694 






2586 65 A 


69.1 


-0.6295 


3.9317 






2600.17 A 


239.0 


-0.0892 


3.5259 


Nin 


-1.4368 


1709.60 A 


32.4 


-2.5756 


-0.0342 






1741.55 A 


42.7 


-2.4477 


-2.4382 






1751.92 A 


27.7 


-2.6331 


-1.2263 


Zn ii 


-2.4081 


2026.14 A 


489.0 


-2.2943 


5.0228 






2062.66 A 


256.0 


-2.5676 


3.2528 



2. A statistical analysis of many-multiplet method 
results 

The central working hypothesis of MWF for evaluating Aa/a 
is the proportionality between the line shift and the atomic rela- 
tivistic correction for all atomic transitions detected in the spec- 
trum of each component of an absorption system. In the frame- 
work of the MM method, one expects the following relation to 
hold for each atomic transition ;: 



Av ; Aa 
— = Qi—+K z 

V; a 



(1) 



where Av,/v, are the line shifts measured for a specific com- 
ponent of the absorption system, with their intrinsic sign. The 
coefficients Q t are given by the ratio qJoJi, where at, is the tran- 
sition wavenu mber, and g, is the a tomic relativistic correction 
(see Table 2 of M urphy et al.ll2003al) . and K z is the redshift of 
that specific component. If Avj/y, and <2, are correlated, MWF 
analysis derives a value of Aa/a different from zero. 

We start with the working hypothesis that there is a pri- 
mary correlation between the relative optical depth of a tran- 
sition and Avj/vj, i.e. that the primary correlation is not that 
between Q, and Av,/v,. A preliminary step is to estimate, for 
a "typical" absorber, the relative optical depths of all transi- 
tions used by MWF. We define the relative optical depth of 
the z'-th transition, r,, as the optical depth scaled to that of a 
reference transition (from now on we skip, for simplicity, the 
index z). In this paper we will use Fe n at 2382A a s the refer- 
ence transition. We take a samp le of low-z absorbers ( Churchill, 
1 19971 IChurchill & VogtL l200ll) to determine the column den- 
sity of Mg i and Mg n relati ve to that of Fe n , and a sample of 
high-z absorbers (Prochaska & Wolfe, 1999; Prochaska et al., 



to determine the column density of Aln, Aim, Sin, 
Cr ii , Ni ii , and Zn n relative to that of Fe ii . We have linked the 
low and high-z sets of data to derive the relative optical depths 
of all transitions of interest. Table[2]gives the decimal logarithm 
of the relative optical depth, log t, the logarithm of the relative 
column density, log iV re i, the Q coefficient, and some other ba- 
sic data ( ionic species, waveleng th, and oscillator strength /, as 
given bv lProchaska et all ( 1200 1) for all transitions used. 

We define the standard deviation in logr, and the linear 
correlation between log r and Q as: 



2 (log T - <l0g T> obs ) 2 



A = 



\ «obs - 1 

E(e-<6)obs)(logT-<logr) obs ) 

(Q - <6>ob S ) 2 V2 Gog T - < lQ g T >obs) 2 ' 



(2) 
(3) 



where the sums are performed only over the n b s transitions in 
which a given absorption system has been observed, and that 
have bee n used to derive Aa/a for that absorption system (see 
Table 3 in Murp hy et al .,2003a, for the sets of transitions used). 
Also the averaged values (logr) b s and (Q) b s are computed 
using only the « D b s transitions that have been observed. 

The values of s and A may change from one absorption sys- 
tem to another, but there are absorption systems with equal val- 
ues of s and A. In principle, spectral lines should be weighted 
with their specific relevance to the fit; however, in our simpli- 
fied analysis we assume that all transitions used for an absorp- 
tion system are equally relevant. The use of log t, instead of 
t, has the advantage of limiting the dispersion, even though the 
resulting correlations are qualitatively similar to those obtained 
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0.2 0.4 0.6 0.8 1.0 1.2 
s = spread in log(r) 

Fig. 1. (a) Distribution of the data in the s-z parameter plane 
and (b) in the s-A plane. Crosses are for objects in the low-z 
sample and circles for the high-z sample (filled circles indicate 
objects in the "high-contrast" subsample). In (b) the number of 
points is smaller than in (a ), because in (b) a single mark may 
indicate more than one absorption system. 



purely based on the set of transitions used, effectively sepa- 
rates lower and higher-z absorbers, with a shallow cut around 
z ~ 1.8. Values of s for low-z absorbers are limited to the range 
(0.3,0.7); values of s for high-z objects lie in the range (0.3,1.2). 
gh-contrast absorption systems typically correspond to higher 
values of s (say, above 0.8): this indicates that, at least on av- 
erage, the quantity s accounts reasonably well for the spread in 
relative optical depth of the atomic transitions used. 

Because different sets of transitions are used for differ- 
ent redshifts, any cosmological implication of measured vari- 
ations of Aa/a with redshift should be taken very cautiously. 
Noti ce that, in a ddition, the low-z and high-z samples given 



by Murphv et a 1. (2003a) correspond to different kinds of ab- 



* ' ' ^ J2003ai) 
sorbers. The low-z sample is mostly composed of absorption- 
line systems selected by the Mg n lines, and have moderate HI 
column densities; while a large fraction of the high-z sample 
is made up of damped Lya absorption systems. Therefore, we 
should not be surprised if, in our analysis, different results are 
found for the two sets of systems, since they have been ana- 
lyzed in different sets of transitions, and furthermore are asso- 
ciated with different types of cosmic structures. 

In the two subsections below we will search for whether in- 
dividual values of Aa/a depend on A. A significant correlation 
would prove that our working hypothesis is correct; if instead 
the MWF results correspond to a true variation of a no correla- 
tion should be found. For A positive, the (spurious) correlation 
between Av/v and Q will have the same sign as the primary 
correlation between Av/v and log r. For A negative the correla- 
tion will have the opposite sign. In order to obtain a negative 
Aa/a when A is positive, Av/v must decrease for increasing 
logr. When A is negative Av/v must increase for increasing 
log r. According to our working hypothesis, the closer the A 
value is to +1, the larger the displacement from zero of the 
Aa/a value for each system. Therefore we expect a correla- 
tion between Aa/a and A. The average value of Aa/a over the 
whole sample will be different from zero, provided that the A 
distribution is not symmetric around zero. 



using r. If log r correlates with Q, the primary correlation be- 
tween log r and the measured Av/v induces a spurious correla- 
tion between Q and Av/v. In the MWF analysis this leads to a 
Aa/a different from zero. 

Our complete s ample is that listed in Table 3 of 
Murp hv et alJ (2003a). We subdivide this into low-z and high- 
z subsamples, using a spectroscopic criterion: an absorber is 
included in the high-z sample only if it has been observed 
in at least one transi tion at wavelengths shorter than 2300A. 
iMurphv et alJll2003ah have noticed that a subsample of systems 
(called the "high-contrast sample"), observed in many transi- 
tions with very different relative optical depths, has a statistical 
spread in Aa/a much larger than the average nominal uncer- 
tainty of each absorber. The authors ascribe this discrepancy to 
uncertainties not accounted for in their analysis, but they ex- 
clude that this may lead to systematic effects in the determina- 
tion of Aa/a. We shall reach a different conclusion. 

Fig- HI shows the complete set of data in the s-z plane, and 
in the s-A plane. Fig. ^a) shows that our criterion, although 



2.1. Low-z sample 

Fig.^7?) shows that negative values of A are mostly associated 
with low-z objects (64 objects). A weighted linear regression 
between Aa/a and A results in a positive slope, although at a 
low significance level: m = (0.80 ± 0.60) x 10 , with a reduced 
X 2 of 1.10 (assuming a functional dependence Aa/a = mA + 
q). A tighter result (4.5 <x) is obtained assuming a functional 
dependence Aa/a = m$A: m = (0.92 ± 0.20) x 10~ 5 , with 
a reduced^- 2 of 1.09. Below we again will use the symbols m 
and mo to indicate the best fitting slope leaving q free, or setting 
q = 0. 

As can be seen in Fig. ^b), low-z objects present a bi- 
modal distribution, consisting of one component with A in the 
range (-0.92,-0.80), and of another one with A in the range 
(-0.52,0.01). The former component is highly uniform in A 
(39 objects with an average A value of -0.85) and the aver- 
age Aa/a is (-0.80 + 0.18) x 10~ 5 . For the other component 
(25 objects with an average A value -0.35) we have instead 
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Aa/a = (-0.36 + 0.23) x 10~ 5 . A comparison of the two sub- 
samples again suggests that Aa/a increases with A. However, 
as already obtained from the linear regression, the significance 
of this trend becomes compelling only if we assume that Aa/a 
vanishes when A = 0. 

It is worth in vestigating another recently published work 
faiand et all 12004). which analyses the spectra of a sample 
of low-z absorption systems. The authors do not find any ev- 
idence of a VFSC and give an upper limit of 0.06 x 10~ 5 to 
Aa/a. This is about one order of magn itude smaller than t he 
signal claimed by M WF. The sample of | Chand et alJ (|2004) is 
smaller than that of Murnh v et alJ (E003a) but the quoted er- 
rors of Aa/a are smaller because high resolution observations 
allowed them to select only systems with all individual com- 
ponents identified unambiguously. The decision of applying a 
severe selection of the systems could either render the approach 
cleaner than that of MWF, or potentially more subject to biases. 
Here we check whether the correlation between Aa/a and A, 
found in the MWF low-z sample, is also present. 

Values of A in Chan d et alJ §004) were obtained using 
all transitions listed for each of the 23 absorption systems 
in their Tabl e 3 (ev en if some components were analyzed by 
IChand et alJ |2004) using a subsample of these transitions). 
The estimated slope for the A- Aa/a linear regression is m = 
(0.15 ± 0.14) x 10~ 5 , with a reduced^ 2 of 0.94; in the case of a 
pure linear proportionality we have mo = (0.12 ± 0.08) x 10~ 5 , 
with a reduced^ 2 of 0.91. This result is consistent with no 
correlation between A and Aa/a, down to a 1.1-1.5 cr level, 
name ly at a level 5-8 times smaller than the correlation found 
in the Murphv *et alJ J200 3a) results. We have also used differ- 
ent methods of object selection, such as limiting our choice to 
systems in which both Fe and Mg transitions are detected in all 
components (12 objects) or to systems in which the same set of 
transitions has been used for all components (8 objects). The 
conclusion, to a rather high l evel of confidence, is that the cor- 
relation is not visible in the IChand et alJ J2004h results. Their 
method does not se em to be affected by th e unwanted system- 
atic effects found in Murp hy et alJ y003a). 



2.2. High-z sample 

Fig-ID^,) shows that high-z absorbers are distributed more con- 
tinuously over a wider range of A, namely (-0.6, +1.0). With 
respect to the low-z sample, the high-z sample shows a much 
larger dispersion in s (see Fig.Q, and therefore it may be better 
suited for searching for effects which arise when using transi- 
tions with a wide range of relative optical depths (i.e. cases with 
large s values). Fig.[5]shows how the measured values of Aa/a 
depend on s in this sample. While no apparent trend is present 
for s < 0.8, for larger s values there is a clear trend of decreas- 
ing Aa/a with s. A difference between these two subsamples 
is present also in the measured values of Aa/a: the average of 
Aa/a for s < 0.8 (38 objects) is (-0.04 ± 0.24) x 10~ 5 , for 
s > 0.8 (26 objects) is (-0.82 ± 0.19) x 10~ 5 . The discrepancy 
between these two values is at a 2.6 cr level. A result similar to 
the s > 0.8 sample, namely Aa/a = (-0.68 ± 0.20) x 10~ 5 , is 
obtained for the high-contrast sample: this is a rather obvious 
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Fig. 2. Dependence on the spread in the relative optical depths 
s of the Aa/a values, for the high-z sample. The dot-dashed 
line divides the figure in two parts: for s < 0.8 no trend of 
Aa/a with s is apparent; for s > 0.8, instead, values of Aa/a 
decrease for increasing s. As in Fig.[T] filled circles indicate the 
"high-contrast" subsample. 



result, since the two samples are partially overlapping. In the 
remaining part of this section, we concentrate on a statistical 
analysis of the high-contrast subsample. 

The high- c ontrast sample is made of 22 objects for which 
iMurphv et alJ (2003a) have measured an excess dispersion of 
Aa/a. We propose that this excess dispersion results from the 
presence of hidden variables, one of which is A. A simple way 
to derive how Aa/a varies with A is to evaluate the average 
Aa/a for systems with negative and positive A: we obtain Aa/a 
equal to (+0.63 + 0.43) x 10" 5 and (- 1.04 ±0.22) x 10~ 5 respec- 
tively, two values which are 3.5 cr apart. A weighted linear re- 
gression between Aa/a and A on the high-contrast sample gives 
a slope m = (-3.46 + 0.88) x 10~ 5 (3.9 cr level). Assuming that 
Aa/a vanishes at A — 0, the slope is mo - (-3.21 +0.62)x 10" 5 
(5.2 cr level). 

In spite of the good significance level of Aa/a- A regres- 
sions for the high-contrast sample, their reduced x 1 values 
are still large , of the order of 3. A possibility is that the 
Murp hy et all <2003al) hypothesis applies mostly to heavily 
damped systems, because they are complex systems, contain- 
ing many saturated lines, while our conjecture of the presence 
of hidden variables may better apply to absorption systems be- 
low a critical column density, N c . 

In Table l2.2l we list all high-contrast systems, and we indi- 
cate the Hi column density whenever this can be found in the 
literature (the relative references are listed in the last column). 
For a few cases the H i column density is unknown: we assume 
for those that the Hi column density is below N c . In Fig.|3jfl) 
we show the data of Table 2 for iV(H i ) > N c = 3 x 10 20 cm" 2 . 
A weighted linear regression for these systems gives a slope 
m = (-1.08 ± 1.15) X 10~ 5 (with a reduced^ 2 of 3.62); or alter- 
natively m = (-1.42 + 0.97)xl0~ 5 (with a reduced^ 2 of3.37). 
Therefore there is no evidence of correlation in this subsample. 
For Af(Hi) < JV C = 3 x 10 20 cirT 2 instead a weighted linear 
regression gives a slope m = (-9.08 ± 1.85) x 10~ 5 (with a re- 
duced^- 2 of 0.91); Fig.|3l^) displays this result. Alternatively, 
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Table 2. High-contrast absorption line systems. 



QSO 




Aa;/a;(10~ 5 ) 


log N(H i ) 


Ref. fl 


0100+13 


2.310 


-3.941+1.368 


21.40 


(PW99) 


0149+33 


2.140 


— 5.1 12zt2.1 18 


20.50 


(PW99) 


0201+37 


1.955 


1.989± 1.048 


20.20 


(PW96) 


0201+37 


2.462 


0.572±1.719 


20.40 


(PW96) 


0528-25 


2.141 


-0.853±0.880 


20.70 


(L96) 


0528-25 


2.811 


0.850±0.846 


21.20 


(L96) 


0841+12 


2.374 


2.277±3.816 


20.95 


(PW99) 


0841+12 


2.374 


1.435± 1.227 


20.95 


(PW99) 


0841+12 


2.476 


-4.304+1.944 


20.78 


(PW99) 


1011+43 


2.959 


2.475± 1.706 




(C02) 


1215+33 


1.999 


5.648+3.131 


20.95 


(PW99) 


1759+75 


2.625 


-0.750± 1.387 


20.80 


(P01) 


1759+75 


2.910 


-0.492± 1.645 


19.80 


(099) 


1850+40 


1.990 


-1.663±0.859 




(C02) 


2206-20 


1.920 


1.878±0.702 


20.65 


(PW97) 


2230+02 


1.859 


-5.407± 1.179 


<20.85 


(PW99) 


2230+02 


1.864 


-0.998±0.492 


20.85 


(PW99) 


2231-00 


2.065 


-2.604± 1.015 


20.56 


(PW99) 


2231-00 


2.065 


1.707± 1.249 


20.56 


(PW99) 


2231-00 


2.653 


-3.348± 1.904 


<20.30 


(L91) 


2343+12 


2.430 


-1.224±0.389 


20.34 


(L98) 


2359-02 


2.095 


-0.068+0.722 


20.70 


(PW99) 



" CL911 lLanzetta et alJ Jl99lh . ( L961 iLu et all <1996t). CL98) 
iLu. Sargent & Barlow! il998h. COggilOutram. Chaffee & Carswelll 
Jl999f). (PW96) IProchaska & Wolfd <1996l). ( PW97 ) 
IProchaska & Wolfd Jl997ft fPW99 ^ IProchaska & Wolfd Jl999h . 
CPOn lProchaska et alJi20oTh . fC02l lCurrant et alJfeCXm 
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m Q = (-4.46 ± 0.81) x 10 -5 (with a reduced X 2 of 1.87). Here 
both regressions determine the slope at a significance level of 
about 5 cr. The result we find is striking: the correlation be- 
tween ha I a and A is very tight for low column density systems 
and it loses its significance for heavily damped systems. One 
worry about choosing N c as low as 3 x 10 20 cm~ 2 might be the 
3 unknown H i column densities of Table 2. These can turn out 
to be higher than 3 x 10 20 cirr 2 . However, we have checked 
that as soon as N c < 7 x 10 20 crrT 2 the significance of the cor- 
relation stays above 4.5 cr for A^(Hi) < N c (even though the 
dispersion increases) and stays below 1 cr for N(H i ) > N c . For 
N c = 7 x 10 20 crrT 2 for example, our last two statements are 
true not only if all the unknown H i column densities are be- 
low N c but also if we do not consider for the correlations the 
two entries relative to the (C02) reference of Table 2, or if we 
consider one or both of them to be higher than N c . 

3. Differential gravitational redshift in absorption 
systems 

In this section we discuss GRS in absorption systems, and its 
possible influence on the quoted variation of the variable fine 
structure constant. GRS can be detected in many types of bound 
systems such as stars, galaxies and clusters of galaxies (e.g. 
iBroadhurst & Scannapiecol 1200(1 IStiavelli & Settilll993b . and 
has been used to infer physical parameters like the total mass. 
It is likely that absorption systems are localized in dark matter 



Fig. 3. Correlations between Aa/a and A for high-contrast 
(high-z) objects, (a) N(Hi) > 3 x 10 20 cirT 2 , (b) N(Ui) < 
3 x 10 20 ctrT 2 . 



halos of galaxies, but since they can be observed only along 
one line of sight, whose position with respect to the system 
barycentre is unknown, it is not easy to quantify the magnitude 
of GRS in observable quantities. To check explicitly if and how 
GRS plays a role in the MWF analysis is even more difficult 
because their data and procedure are not available to the public. 

In order for GRS to be responsible for the fine structure 
constant variations measured in absorption line systems with 
the MM method, three conditions are required: 

i) the differential GRS between the innermost regions and the 
outskirts of typical absorption systems should be at least of the 
same order as the relative shifts measured by MWF (of the or- 
der of 0.1 km s , and interpreted in terms of a VFSC); 

ii) misidentification of individual components in a line profile 
should affect the MM analysis, as carried on by MWF, i.e. com- 
ponents attributed to the same cloud in different transitions are 
instead associated with different clouds; 

Hi) the kinematics should allow a detection of a tiny spectral 
shift of gravitational origin. 

In the rest of this Section we examine these conditions in 
more detail. In Section 4, we test in a more quantitative way 
the GRS hypothesis, by searching for correlations between the 
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estimated Aa/a and indicators of the gravitational potential in 
the absorption system. 

3.1. Galactic haloes and gravitational redshift 

The aim of this subsection is to evaluate, for typical gaseous 
systems embedded in dark matter haloes, the magnitude of the 
differential GRS. Photons of frequency vo emitted in a gravita- 
tional potential <I> are observed at infinity to be redshifted by an 
amount 



Av 



75' 



(4) 



where Av = v - vo, and c is the speed of light. 

Let us consider gaseous systems embedded in dark matter 
haloes whose density profil e is described by th e formula pro- 
nosed bvlNavarro. Frenk & Whird(ll99fiiri997l): 



p(R) = 



Pnfw 



(R/R NFV /)(l + R/R NFW ) 



(5) 



If R& is the radius of a sphere containing a mean den- 
sity A times the critical density at redshift z, and Ma is the 
mass inside R&, one can express the gravitational potential 
of a dark matter halo at a redshift z in terms of Ma and of 
the c o ncentr ation parameter Ca = /?a/^nfw JNavarro et all 
1 19961 Il997h . Numerical simulations of structure formation 
often reveal a correlation between Ca and Ma which de- 
pends on the details of die assumed cosmological model (e.^ . 



Navarro, Frenk & White, 1996, 1997; Avila-Reese et 



toael (e.g. 

■ail boou 



Eke, Navarro & Steinmetz, 2001). For a Cold Dark Matter- 



dominated universe, the dependence of Ca on Ma is weak. 
We consider A = 200, as originally proposed bv lNavarro et alJ 
ill996Ml997l) . and define C = C2oo- Here we do not assume 
any explicit cosmological model to relate C and M200 but use 
C ^2-10 at redshifts 0.5-3, as suggested by numerical simu- 
lations of dark matter haloes forming in a Cold Dark Matter- 
dominated universe. The gravitationally redshifted frequency 
of an absorption line inside a galactic potential dominated by 
dark matter can therefore be written as: 



Av 



GM200 ln(l+Cx) 
~4jtc**200 41n(l + C) - C/(l + C)} 
2.3 x 10- 8 (M 2 oo/h- l \ 213 (1 +z)ln(l + Cx) 



(Q Q) 1/3 \ 10 12 M Q J x{ln(l + C) - C/(l + C)} 



(6) 



where G is the gravitational constant, Q is the matter den- 
sity of the Universe at z and we have introduced the variable 
x = R/R20Q. In Fig.|4]we show the GRS in km s _1 as a function 
of the galactocentric distance x at z = 2 for Qo/Q - 1. For a 
flat, cosmological constant dominated universe with £2o = 0.3, 
Ao = 0.7, the GRS is a factor ~ 1.5 larger than what Fig.|4] 
shows. There are no noticeable differences in the magnitude 
of the GRS if we consider other gravitational potentials with 
simil ar masses but no cusp in the centre. For a Burkert poten- 
tial jBurkertL 1 19951) . for example, we find less than a factor 2 
difference over the whole range of masses and radii shown in 

Kg. a 

An interesting result of Fig.|4]is that the velocity shifts are 
of the same order of magnitude as those found in the analysis 
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Fig. 4. GRS in kms 1 as a function of the galactocentric dis- 
tance x = R/R200 at z - 2. The 3 sets of curves correspond to 
the halo masses M 20 oh = 10 n , 10 I2 ,2 x 10 12 M . The contin- 
uous lines are for concentration values C = 2, the dashed lines 
for C = 6, the dash-dotted lines for C = 10. 



of Murph y et all d2003al see their Table 4) and used to claim a 
VFSC. The estimate of differential GRS presented above refers 
to the whole absorption system; differential GRS within indi- 
vidual clouds (detected as individual components in the line 
profile) is much smaller. 

3.2. Detectability of GRS in QSO absorption systems 

We outline here the most appropriate conditions to detect GRS 
effects in QSO absorption systems. In principle, GRS could be 
detected using even a single transition, if the available sample 
of absorption lines is large enough to justify a statistical ap- 
proach. In fact, isotropy arguments show that, in the absence 
of GRS, kinematic redshifts of individual spectral components 
should average to zero with a symmetric distribution in the 
whole sample. GRS, instead, induces skewness in the spectral 
distribution of the absorption components, because absorption 
lines located more deeply in the potential well are those more 
strongly redshifted. In this kind of analysis, the main practi- 
cal difficulty is to cope with the smallness of this effect, com- 
pared to both instrumental and statistical uncertainties. Even in 
the absence of instrumental uncertainties, there is a minimum 
number of spectral components below which it is impossible 
to extract the gravitational shift from the statistical dispersion 
due to kinematic motion. The problem is qualitatively similar 
to that of determining, with accuracy AVV, the average shift 
for a sample of components whose distribution shows a disper- 
sion AV ± : in this case at least ~ (AV ± / AVo-) 2 different compo- 
nents are required. If we aim at detecting AVV ~ O.lkms -1 , 
while the dispersion of the individual absorption components 
is AV ± ~ lOOkms -1 , at least ~ 10 6 components are required, 
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Fig. 5. Relative change of V glav and V or b, r , along the line of 
sight, for different values of the impact parameter d. 



a number far beyond that available in the present data sets. 
This argument, even though in a more complex formulation, 
should apply also in a multi-line approach. A smaller number 
of components may be required if the velocity pattern of the 
material changes smoothly with position. For illustration, we 
present here a reasoning based on the assumption that all blobs 
move in circular orbits, seen edge-on, in a potential consistent 
with the profile given in Eq.[5] In this case, at a given radius R, 
the orbital velocity and the gravitational redshift (expressed as 
a velocity) are given respectively by: 

y orb = 47r Gp NF wflNFW 

(ln(l + R/R NFW )/R - 1/ORnfw + R)) , (7) 
V grav = 4nGp NFW Rl FW ln(l + R/R NFW )/(Rc). (8) 

Let us also assume that, at R = /?nfw, the ratio V grav /V or b = 
1.577 -y/47rGpNFW^NFW / c is of the order of the observed ratio 
AV, r /AV ± ~ 1(T 3 . 

We consider a path along a line of sight with impact param- 
eter d with respect to the centre of the system: R = V/ 2 + d 2 , 
where I is the coordinate along the line of sight. We then com- 
pute Vgrav an d the radial component of V or b ( V rb,r)- Results are 
displayed in Fig. [5] for different values of d. While by a pure 
dimensional argument one would always expect values of the 
order of AVo-/AV ± ~ 10 3 , Fig.|5]shows that there may be con- 
ditions in which this ratio can be one, or even two orders of 
magnitude higher. In addition, in the case of orbits tilted by 
an angle f , a further 1 / cos(^) factor should be included. This 
shows that, under favorable conditions, a sample of ~ 100 com- 
ponents, spread over a set of absorption line systems, may be 
sufficient to get a statistically significant detection of GRS. 



3.3. Component misidentification: how GRS breaks 
the symmetry 

MFW aim at comparing components, in different transitions, 
which are attributed to the same cloud, i.e. to the same physi- 
cal region within the absorption system. Of course, if this goal 
is reached, no GRS would be detected. Let us examine instead 
the case of "component misidentification", namely the case in 
which two components, identified in the line profile of two 
or more different transitions, are erroneously attributed to the 
same cloud. There may be different reasons why a given tran- 
sition effectively samples inner region clouds compared to an- 
other transition. This may happen, for instance, if chemical or 
ionization radial gradients are present in the absorption system. 
There is some evidence of chemical and ionization inhomo- 
geneities in high column density absorbe rs, for example, w hich 
are very numerous in our high-z sample l|Petitiean & SrianandL 
ll999HPetitiean. Srianan d & Cedricl l2002l) . But even in the ab- 
sence of such gradients, inner regions are usually denser than 
outer regions, and absorption lines of different relative optical 
depth saturate at different distances from the centre and lines. 
For this reason, lines with different t might effectively sam- 
ple different regions in the potential well, and therefore may 
be subject to different levels of gravitational redshift. From 
this scenario it naturally follows that one should find a correla- 
tion between Ay/ v and t for different transitions, and therefore 
a correlation between Aa/a and A (see Section 2). Therefore 
our working hypothesis of Section 2 (a primary correlation be- 
tween Av/ v and r) is derived naturally from the idea that GRS 
and component misidentification may play a role in some deter- 
mination of Aa/a. In our framework, if transitions of different 
relative optical depths sample different regions of an absorp- 
tion system, high-contrast absorption systems should in fact be 
the most appropriate ones to reveal effects of differential grav- 
itational redshift. The correlations we discuss in Section 2.2 
indicate indeed that strong systematic effects are present in this 
subsample. 

In the presence of component misidentification, if clouds 
seen in absorption in two different transitions are typically lo- 
cated at different distances from the centre, their frequency ra- 
tio is affected by the different radial velocities as well as by 
the different halo gravitational potential felt by the two clouds. 
Symmetry arguments imply that spurious deviations caused by 
the different kinematic properties of the clouds average to zero 
in a suitably large sample. Therefore only GRS can reasonably 
account for any residual systematic shift in the average value 
of Aa/a. 

4. Correlations between line frequency shifts and 
gravitational potential 

In Section 2 we have shown that the concept of A is very ef- 
fective in searching for hidden correlations. This does not nec- 
essarily imply that GRS plays a role in generating systematic 
effects, even though some working hypotheses suggest that it 
might. In this section we assume that, for a uniform sample 
in A, the derived value of Aa/a is roughly proportional to the 
level of GRS. Since GRS is a function of the galactic potential, 
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we investigate whether Aa/a correlates with indicators of the 
gravitational potential strength. The gravitational potential de- 
pends on the total mass associated with the absorption system 
and on the radial distance of the absorbing material from the 
centre of mass, but unfortunately none of these physical quan- 
tities is directly measurable from absorption data. A way out 
is to identify observable quantities, which correlate with total 
mass and radius. 

A quantity that can correlate with the strength of the gravi- 
tational potential is the metal column density. In fact, for con- 
stant metal abundances, more massive absorption systems have 
higher gas column density (and therefore higher metal column 
densities). In addition, metals are expected to be overabundant 
towards the centre of the potential well or in more massive sys- 
tems because of the longer star formation history. 

In this section we examine two samples of objects, both 
with known Fen column densities (hereafter A^(Fen)). The 
fir st one contains low-z a bsorption systems (as listed in Table 3 
of IMurphv et all l2003al) that have well documented N(Fen) 
along the line of s i ght. M ost of thes e objects have been ana- 
lyzed bv lChurchilllJl997l C97), or bv lChurchill & Vogtlll200l[ 
CV01). In order to obtain a uniform sample, we have limited 
our selection to these two references. The resulting sample 
contains 24 absorbing systems at redshifts 0.5 < z a bs < 1-8, 
and their logA^(Fen) is listed in Table 3. This sample is 
statistically consistent with the larger low-z sample analyzed 
in the previous section. Our second sample contains high-z, 
damped Lya absorption systems (with A^(H i ) > 10 20 cirT 2 and 
Zabs < Zem)- In order to have a nearly uniform sample we select 
22 of them with 1.8 < z a bs < 3.5, and published values of 
log jV(Fe ii ) in at least o ne of the following referenc es: (L96) 
Lu Let all liooS. (PW96^ IProchaska & Wolfd JT996T). (PW97^ 
Prochaska & Wolfe J dl997l). (PW9 9) IProchaska & W olfe 
\ 1 999l) . (P0 1 ) IProchaska et al.1 d200lh . The sample is listed in 
the last 21 entries of Table 3. 

In order to estimate the dependence of the measured values 
of Aa/a on A^(Fen), we will use a linear regression, adopting 
the functional dependence Aa/a = ra(log A^(Fe n ) - 13.5) + q. 
In this way, for the whole low-z sample of Table 3, we obtain 
m = (-0.45 ± 0.33) x 10~ 5 and q = (-0.50 ± 0.25) x 10~ 5 , 
with a reduced x 1 of 0.90: the value obtained for the the slope 
m is not statistically significant (1.4 cr). For the high-z sample 
of Table El instead, we calculate m = (1.19 ± 0.51) x 10~ 5 and 
q = (-2.21 ± 0.72) x 10~ 5 , with a reduced^- 2 of 2.65: a rea- 
sonable (2.3 cr) significance level is obtained here, but the re- 
duced^- 2 is anomalously high. We attribute this result to the fact 
that some of the quoted uncertainties on Aa/a, which are rela- 
tive to high-contrast damped systems, are underestimated (see 
Fig.fJJaJ, and the discussion in the previous section). In fact, if 
we compute a linear regression on the high-z sample, assigning 
the same uncertainty to all measurements of Aa/a, we obtain 
m = (-0.38 + 1.20) x 10~ 5 and q = (-0.27 ± 1.45) x 10~ 5 : 
therefore, also in this case the slope m is not statistically signif- 
icant. We attribute the low significance of these correlations to 
the fact that none of the two samples is uniform in A. 

Luckily, 18 out of the 24 objects of the low-z sample lie 
in a narrow range of A, namely (-0.92, -0.79). The correlation 
becomes more significant if we select only absorbers with high 



Table 3. Fe n -selected absorption line systems. 



QSO 




Ao-MlO" 5 ) 


logjV(Fen) 

CUT 2 


Ref." 


0002+05 


0.851 


-0.346± 1.279 


13.802 


(CV01) 


0117+21 


0.729 


0.084+1.297 


12.839 


(C97) 


0117+21 


1.048 


-0.223±2.200 


12.718 


(CV01) 


0117+21 


1.343 


-1.290+0.948 


12.484 


(C97) 


0420-01 


0.633 


4.211±4.076 


12.879 


(CV01) 


0450-13 


1.175 


-3.070±1.098 


14.985 


(C97) 


0454+03 


0.860 


0.405± 1.325 


15.068 


(CV01) 


0454+03 


1.153 


-0.749+1.782 


12.826 


(CV01) 


0823-22 


0.911 


-0.394±0.609 


13.597 


(CV01) 


1148+38 


0.553 


-1.861+1.716 


13.186 


(CV01) 


1206+45 


0.928 


-0.218+1.389 


12.905 


(CV01) 


1213-00 


1.320 


-0.738+0.760 


14.529 


(C97) 


1213-00 


1.554 


-1.268+0.892 


14.404 


(C97) 


1222+22 


0.668 


0.067+1.474 


13.115 


(CV01) 


1225+31 


1.795 


-1.296+1.049 


13.881 


(C97) 


1248+40 


0.773 


2.165+1.191 


13.550 


(CV01) 


1248+40 


0.855 


-0.021+1.268 


12.533 


(C97) 


1254+04 


0.519 


-3.371+3.247 


13.957 


(CV01) 


1254+04 


0.934 


1.485+1.908 


12.699 


(CV01) 


1317+27 


0.660 


0.590+1.515 


13.100 


(CV01) 


1421+33 


0.843 


0.099+0.847 


13.143 


(C97) 


1421+33 


0.903 


-0.998+1.783 


13.706 


(CV01) 


1421+33 


1.173 


-2.844+1.448 


13.049 


(CV01) 


1634+70 


0.990 


1.094+2.459 


13.059 


(CV01) 


0000-26 


3.390 


-7.666+3.231 


15.146 


(PW99) 


0019-15 


3.439 


0.925+3.958 


14.770 


(PW99) 


0100+13 


2.310 


-3.941+1.368 


15.096 


(PW99) 


0149+33 


2.140 


-5.112+2.118 


14.202 


(PW99) 


0201+37 


2.462 


0.572+1.719 


15.060 


(PW96) 


0216+08 


1.768 


0.044+1.235 


14.530 


(L96) 


0347-38 


3.025 


-2.795+3.429 


14.503 


(POD 


0528-25 


2.141 


-0.853+0.880 


14.940 


(L96) 


0741+47 


3.017 


0.794+1.796 


14.041 


(P01) 


0841+12 


2.476 


^.304+1.944 


14.434 


(PW99) 


1215+33 


1.999 


5.648+3.131 


14.648 


(PW99) 


1425+60 


2.827 


0.433+0.827 


14.480 


(L96) 


1759+75 


2.625 


-0.750+1.387 


15.091 


(P01) 


1946+76 


2.843 


-4.743+1.289 


13.380 


(L96) 


2206-20 


1.920 


1.878+0.702 


15.458 


(PW97) 


2206-20 


2.076 


1.429+3.022 


13.320 


(PW97) 


2230+02 


1.864 


-0.998+0.492 


15.188 


(PW99) 


2231-00 


2.066 


-2.604+1.015 


14.750 


(PW99) 


2348-14 


2.279 


1.346+4.180 


13.792 


(PW99) 


2359-02 


2.095 


-0.068+0.722 


14.507 


(PW99) 


2359-02 


2.154 


4.346+3.338 


13.895 


(PW99) 



" (C97) Churchill (19971) . (CV0 1) IChurchill & Vod | 200l|) . 
(L96) LuetalJ <1996h. (PW96) IProchaska & Wolfd Jl996l) 
(PW97 ) IProchaska & Wolfd j 19971) (PW99) IProchaska & Wolfd 
Il999l) . (POll lProchaska et alJfcOOlh 



logAf(Fen) values in this A range. Taking logA^(Fen) > 13.0 
(12 objects) we obtain m = (-1.23 ± 0.56) x 10~ 5 (2.2 cr) and 
q = (-0.35 + 0.34) x 10~ 5 , with a reduced x 2 of 1.14. For 
logAT(Fen) > 13.5 (6 objects) m = (-2.16+0.78)xl0~ 5 (2.8 cr) 
and q = (0.34 ± 0.53) x 10~ 5 , with a reduced X 2 of 0.92. This 
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result is consistent with a Aa/a that deviates significantly from 
zero only for logiV(Feii) larger than about 13.5. The Aa/a 
deviation is towards more negative values as logiV(Fen) in- 
creases. Although the statistical significance is not high (only 
2.8cr in the best case), one should notice that negative devia- 
tions of Aa/a for negative A are consistent with the increase 
of Aa/a with increasing A, obtained in the previous section for 
the low-z sample. This can be explained in terms of an effec- 
tive line segregation, with weaker lines sampling effectively re- 
gions closer to the centre of the potential well. GRS is a reason- 
able candidate to explain the decrease in Aa/a with increasing 
log iV(Fe ii ), because systems with higher metal column den- 
sity are either closer to the centre of the potential well, or lo- 
cated in more massive systems or both. 

There have been other proposals to explain the VFSC re- 
sult. Of particular interest is the consistency of the VFSC val- 
ues with a non-solar isotopic ratio of ( 25 Mg + 25 Mg)/ 24 Mg oc- 
curring at large redshifts. If this isotopic ratio increases at large 
redshifts, small apparent shif ts would be introduced in the ab- 
sorption lines ( Ashenfelter et al., 2004). These shifts mimic a 
VFSC if instead a solar isotopic ratio of magnesium is used, as 
in the MWF analysis. The large uncertainty about this proposal 
lies in the magnesium isoto pic ratio variations with metallic- 
itv. lAshenfelter et ai1 (E004) show that an Initial Mass Function 
(IMF) particularly rich in intermediate mass stars is needed 
for an increase of the magnesium isotopic ratio with decreas- 
ing metallicity (in th e range [Fe/Hl-0 , —1.5) . Analyses of low 
metallicity star data l lGav & Lam bert. 2000) show instead an 
opposite trend of the ma gnesium isotopic ratio v ariations with 
metallicity. If the data bv ldav & Lambertl lEoOO) reflect the sit- 
uation in QSO absorbers, which have lower metallicity than 
solar, then magnesium isotopic ratio variations cannot explain 
the negative Aa/a reported by MWF. In this case metallicities 
lower than solar would correspond to lower magnesium iso- 
topic ratios and to posi tive value of Aa/a. If o ne assumes in- 
stead that the results of Ashenfelt er et alJ (|2£)04) apply to gas in 
QSO absorbers, then magnesium isotopic ratio variations can 
explain the MWF results. The correlation between Aa/a and 
metallicity should in this case be positive for a certain range of 
[Fe/H] (which depends on IMF and yields). For the GRS hy- 
pothesis instead large variations of Aa/a should be associated 
with large potential wells. Since more massive galaxies have 
had a longer star formation history, a negative correlation be- 
tween metallicity and Aa/a should appear. The results of this 
Section seem to favour the GRS hypothesis, but accurate metal 
abundances, higher significance in the outlined correlations and 
more information on earlier IMF are needed to draw any con- 
clusions on the magnesium isotopic ratio variation hypothesis. 

5. Conclusions 

Even t hough we cann ot disprove completely the VFSC hypoth- 
esis o f Murphv et alJ d2003ah . in this paper we have discovered 
systematic effects which can mimic a non-zero Aa/a. Often 
Aa/a depends on A, the correlation coefficient between two 
atomic quantities: the relative optical depth and the relativistic 
correction coefficient Q. For the atomic transitions used in the 
MWF analysis for each absorber, we evaluate A and notice that 



the A distribution is not symmetric around zero: this justifies 
the net displacement of the average Aa/a from zero. In par- 
ticular, the correlation between spectral shifts and the relative 
optical depth of the various transitions seems to be the origin 
of correlations between spectral shifts and Q, used by MWF as 
evidence for a VFSC. Non-zero values of Aa/a may then be 
the result of GRS if lines of different relative optical depth ef- 
fectively sample regions at different distances from the centre 
of the potential well. 

We have e xamine d the low-z and the high-z samples of 
iMurphv et alJ J2003al) separately, since they have been ana- 
lyzed in different sets of transitions and furthermore are asso- 
ciated with different types of cosmic structures (see Section 2 
for more details). In the low-z sample the correlation between 
Aa/a and A can be interpreted in terms of a primary positive 
correlation between Av/v and the relative optical depth: that 
means that relatively weaker lines are more redshifted than oth- 
ers (i.e. Av more negative). This is in agreement with what one 
would expect from GRS, under the hypothesis that intrinsically 
weaker lines are visible only through dense regions, located 
closer to the centre of the absorption system. Since for this sam- 
ple most of the systems lie in a restricted A range, and we know 
the metal column densities, we compare iV(Fe n ) with the de- 
rived Aa/a values. Systems with higher A^(Fen) values have 
more negative Aa/a values. This again confirms the GRS hy- 
pothesis since more negative Aa/a can be interpreted as more 
heavily gravitationally redshifted lines. High A^(Fen) values 
are in fact expected when the gas lies in regions of strong grav- 
itational potential. If future data prove a positive correlation 
between Aa/a and metal abundances, this will support GRS- 
related studies and weaken other proposals, such as the magne- 
sium isotopic ratio variations, for explaining the non-zero value 
of Aa/a derived in the MWF analysis. 

For th e "high-contrast" subs ample of the high-z systems, 
for which IMurphv et alJ J2003al) have reported an anomalous 
statistical dispersion in Aa/a, we find instead a systematic ef- 
fect that can be interpreted in terms of a primary inverse cor- 
relation between Av/v and the relative optical depth. Further 
dividing this subsample into heavily damped systems, with 
N(Ui) > Af c , and systems with JV(Hi) < A^ c , we have dis- 
covered that the latter subsample shows a striking inverse cor- 
relation, at a level higher than 4.5 cr for N c < 7 x 10 20 cm" 2 . No 
correlation is present for A^(H i) > N c . So, while for the heavily 
damped systems the MWF interpretation of the anomalous sta- 
tistical dispersion might be correct, for lower column density 
systems we believe that the anomalous statistical dispersion is 
indeed the effect of an underlying correlation. 

Note that the sign of the correlation between Aa/a and A 
for the high-contrast (high-z) sample is opposite to that ob- 
tained for the low-z one. If this correlation is interpreted in 
terms of GRS, stronger transitions should be associated with 
innermost regions of the absorption system. This is opposite to 
what is inferred for the low-z sample, and requires the pres- 
ence of physical inhomogeneitie s within the system ( as ac- 
tually found in some cases, e.g. Petitiean & Sria nandl 1 19991: 
IPetitiean et all l2002l) . However, one cannot exclude that the 
opposite trends simply arise because the bias comes out dif- 
ferently in the MWF procedure used to estimate Aa/a: this is 
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because the data and the characteristics of the absorption sys- 
tems in the two samples are very different. 

As already pointed out, a more detailed analysis can be 
carried out having access to the original data. For instance, 
it would be useful to obtain subsamples analyzed uniformly, 
using the same set of transitions. Conversely, one could also 
investigate how Aa/a changes with different choices of the 
fitted transitions, in the same absorption system. More gener- 
ally, instead of determining how Av/ v correlates with Q (from 
which Aa/a is determined), one should search for correlations 
of Av/v with other quantities, such as the relative optical depth. 
Notice that any quoted variation of Aa/a with cosmological 
time should be taken very cautiously, because different sets of 
transitions are used for different redshifts. 

In a recent work, IChand et alJ J20041) found no variation 
in the fine structure constant for a sample of low-z absorption 
systems. They obtained this result avoiding both weak and sat- 
urated lines, as well as heavily blended spectral regions. In this 
way they restricted their analysis mostly to intermediate "satel- 
lite components", avoiding the ones relative to central or outer- 
most regions. Therefore they reached a better accuracy in posi- 
tioning the subcomponents. A contamination of GRS on Aa/a 
is possible only in the presence of some bias in the spectral 
analysis. We do n ot find any signific ant correlation between 
Aa/a and A in the lChand et alJ d2004l) results. Their work fur- 
ther supports the idea that biases may originate from the fit 
of weak spectral components or of the saturated and complex 
parts of a line structure. Moreover, GRS effects are harder to 
detect in a sample which avoids the centremost and outermost 
regions since this choice effectively restricts the range of the 
potential well tested. 

We would like to encourage any future GRS experiment 
over cosmological distances since the results can be unique 
tools to test structure formation scenarios in the context of sev- 
eral cosmological models. For example if at redshifts as high as 
4-5 the Universe is still dominated by small mass objects (Cold 
Dark Matter scenario), GRS-induced line frequency shifts in 
absorption systems at such redshifts should be smaller than the 
ones detectable at more recent times. Also theoretical models 
of possible variations of the fine structure constant, which take 
into account the ob servational state of the art are needed (see 
e.g. lSteinhardl [20Q3>. in order to check the consistency of dif- 
ferent results. 
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